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Abstract 

Dalitz decay of baryon resonances is studied and expressions for the 
decay width are derived for resonances with arbitrary spin and parity. 
Contributions of the various terms in the transition matrix element 
are compared and relevance of spin-parity and the resonance mass 
is discussed. Explicite algebraic expressions are cited for spin<5/2 
resonances. The results can be used in models of dielectron production 
in elementary reactions and heavy ion collisions. 
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1 Introduction 



The aim of relativistic heavy ion collision experiments is to study the prop- 
erties of the hot and dense nuclear matter created in the collisions. These 
properties have to be reconstructed using the information encoded in the 
particles detected in the final stage of the collision. The best probes for this 
investigation are dileptons (in the few GeV per nucleon bombarding energy 
range dielectrons) , which leave the hot and dense zone unaffected by final 
state interactions. 
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Theoretical studies predict a modification of hadron properties (masses 
and decay widths) in hot and dense nuclear matter. These medium effects 
can most easily be studied in the case of vector mesons, which decay directly 
to dileptons and their masses are measured through the invariant mass of 
the resulting dilepton. 

Dielectron production in heavy ion collisions with a few GeV per nucleon 
energy has been studied experimentally by the DLS collaboration at BE- 
VALAC and will be studied by the HADES at GSI. In order to obtain a 
reliable model of heavy ion collisions first one has to check the elementary 
reaction channels used in the model. For this purpose dilepton production 
in nucleon-nucleon collisions has to be investigated both experimentally and 
theoretically. 

Particle production in nucleon-nucleon collisions is usually described us- 
ing the resonance model, i.e., the final state particles are produced via the 
creation and subsequent decay of baryon resonances. This multistep process 
is factorized so that initial and final state particles of the individual steps 
are handled as on-shell particles and in case of resonances their masses are 
generated according to a Breit-Wigner distribution. Different factorizations 
result in different models with different elementary channels. 

In models of dielectron production in pp collisions usually the following 
elementary channels are included: direct dileptonic decay of p°, u and <j) 
vector mesons, Dalitz decay of mesons tt° — > / -/e + e~ , 77 — > r ye + e~ , u — > ir e + e~ 
and Dalitz decay of the A(1232) resonance, A — > Ne + e~ . Such models (e.g. 
[0], 0, [|) can reasonably well reproduce the experimental dilepton invariant 
mass spectra in the Eb = 1 - 2 GeV bombarding energy range Q . In models 
of heavy ion collisions the additional channels of pion annihilation and pn 
bremsstrahlung appear. 

With increasing beam energy creation of higher mass baryon resonances 
becomes more and more relevant. In heavy ion collisions kinetic energy of 
colliding particles can be accumulated in multiple scattering processes, there- 
fore higher mass resonances appear at lower beam energies. These resonances 
may contribute to the dilepton spectrum through their Dalitz decay similarly 
to the A (1232). This contribution may especially be relevant in the higher 
dilepton mass range where the p° and A(1232) channels are already negli- 
gible. This dilepton mass range coincides with the range of the <fi meson 
mass. 

Recently, production of mesons in subthreshold heavy ion collisions has 
received some attention both experimentally || and theoretically || [7]. (f> 
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mesons are being studied by the FOPI group at GSI via their kaon decay 
channel, <fi — > K + K~, but the dileptonic decay channel will also be accessible 
by the HADES detector. In the latter case an important source of background 
may come from the Dalitz decay of higher mass resonances. 

In this paper we derive expressions for the differential Dalitz decay width 
of baryon resonances with arbitrary spin-parities. We discuss in detail the 
Dalitz decay of resonances with spin<5/2. Our paper is organized as follows: 
in Sec. [| we derive expressions for the electromagnetic transition matrix el- 
ements of baryon resonances with arbitrary spin-parity and give expressions 
for the Dalitz decay width. In Sec. we give a detailed discussion of the 
numerical results. We compare contributions of the various terms in the 
transition matrix element and discuss the relevance of spin-parity and the 
resonance mass. In the Appendix we cite the spin projectors used in the cal- 
culations and the explicite expressions for the polarization averaged squared 
transition matrix elements of the baryon resonance Dalitz decays. 

2 Dalitz-decay of a spin- J resonance 

2.1 Matrix elements of the electromagnetic current 

To study the electromagnetic decays of a spin- J baryon resonance one has 
to determine the matrix element of the electromagnetic current operator J^ 
between the nucleon (N) and the spin-J baryon (R) state {N\J^\R). Let p*, 
m*, and A* denote the four-momentum, mass, and helicity of the R resonance, 
respectively, and p, m, and A the corresponding quantities for the nucleon. 
Let us introduce the notations q = p* — p and P = (p* + p)/2. 

For J > 3/2 the spin-J fermion can be described by a generalized Rarita- 
Schwinger spinor-tensor field ^fP 1 "'P n (n = J — 1/2), therefore the momentum- 
space wavefunction of the particle is the spinor-tensor amplitude u Pl "' pn (p*, A*) 
This quantity in general contains some additional lower spin components. 
The spin-J content of u pl ' Pn (p*, A*) can be selected by prescribing the gen- 
eralized Rarita-Schwinger relations: 
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The matrix element of the electromagnetic current can be written generally 

as 

(N\J P \R) = u(p, \)Y lipi ... Pn u pl " pn (p*, A*), (2) 

where the form of T ppi ... Pn is restricted by the conservation of electric charge, 
q fl (N\J p \R) = 0. The relations (0) and the Dirac equations u(p, \){j>—m) = 
and (yf« — m*)u pl '" Pn (p#, A*) = reduce further the number of independent 
terms in T ppi ... Pn to three, which can be chosen of the form 

3 

r MPi-Pn = I] fM 2 )x l Wl Vp2 ■ ■ -P P nG, (3) 

with 

<?M<?p - <?Vp' ( 4 ) 

and G = 1 or 75 for resonances with positive or negative normalities, respec- 
tively. The normality of a spin-J baryon is defined by P(— l) J ~ l l 2 with P 
the intrinsic parity. fi(q 2 ), i — 1, 2, 3 are three independent form factors. 

In the spin-1/2 case the general form of the matrix element of the elec- 
tromagnetic current is 

(N\J P \R) =u(p,\)T IM u(p„K), (5) 
and there are two independent terms in T M : 

r, = X>(g 2 )x;G, (6) 

i=l 

with 

Xi = {P-q)l,-iP^ 

xl = (i 2 ip-i%- (7) 

G is again 1 or 75 depending on the normality of the resonance. 



XfjLp 
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x 3 
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2.2 The Dalitz-decay width 

The differential width of the Dalitz-decay of a particle is related to its pho- 
tonic decay width to a virtual photon, Tr- > n 1 (M), by || 

d^R^Ne+e- Qi 1 _ /ji,r\ fo\ 

dM* = 3^M^ r ^ (M) - (8) 

Here the notation M 2 (= q 2 ) is used for the square of the dilepton invariant 
mass (= mass of the virtual photon). Fr^,^(M) can be expressed in terms 
of the photonic decay matrix element (N^\T\R) as 

J\(mlm 2 ,M 2 ) l ^ _ 

Wmh V _ — giwW. (9) 

where the transition matrix element is related to the electromagnetic current 
matrix element by 

(N>y\T\R) = -e»(N\J,\R), (10) 

with e M the photon polarization vector. In (§) n po i,R is the number of polar- 
ization states of the R resonance and A (a, b, c) = a 2 + b 2 + c 2 — 2(ab + be + ac) 
is the usual kinematical factor. 

The polarization sum in (§) runs over all physical polarization states of the 
incoming and outgoing particles, that is, in the case of the spin- J resonance 
only over those states that belong to the spin- J content of u pi '" Pn (p*, A*). 
As a result, the spin- J projector appears in (0) after substituting (Q). The 
projectors for spin-3/2 and 5/2 are given in Appendix A. 

If the differential cross section of the production of the resonance R in 
an elementary reaction, dan/ dm*, is known Eq. (f$D can be used to obtain 
the cross section of dilepton production through the Dalitz decay of R in the 
same reaction. This differential cross section is given by 

d(7 e +e-,R Damz _ f dm da R 1 dT R ^ Ne+e -(m t ) 

dM J *dm,Tg t (m,) dM 1 ' 

2.3 Form factors and coupling constants 

A possibility to determine the functional form of the form factors fi(M 2 ) 
is the application of the vector dominance model (VDM) ||. However, we 
aim at a model that handles the theoretically more interesting vector meson 
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direct dilepton decays separately from the baryon resonance Dalitz decays 
serving as a background. Applying the VDM to the baryon resonance Dalitz 
decay channels simultaneously with the inclusion of the vector meson direct 
decays would give rise to a double counting if - following the usual treatment 
- the model incorporates vector mesons originating from decays of baryon 
resonances. The authors of [10] use a sophisticated method to avoid this 
double counting. In the present paper we do not apply the VDM and keep 
the form factors fi{q 2 ) constant. 

We introduce the dimensionless coupling constants gi via the equation 
fi = egi/m kt , where e is the elementary charge related to the fine structure 
constant by e 2 = Ana. The exponent ki of the nucleon mass is determined 
by the requirement that the couplings g>j are dimensionless. The values of 
are obtained using the real photonic decay width of the resonances. 



3 Discussion of the results 

We first compare the contributions of the various terms in the matrix ele- 
ments (0) and (^) to the dilepton invariant mass spectrum. For this pur- 
pose we have set the dimensionless couplings to unity and calculated the 
contributions of the individual terms and the interference terms to the dif- 
ferential Dalitz decay width. We carried out these calculations for spin<5/2 
resonances and chose a resonance mass of 1.5 GeV. The resulting plots are 
shown in Fig. [I] 

We notice that the contribution of terms containing g2 in the spin- 1/2 
case and terms containing g^ in the spin>3/2 case is free from singularities 
at M 2 = 0. The real photonic width of the resonance r^_ > jv 7 (M = 0) ob- 
tained from these terms is zero in accordance with Eq. (^), consequently, the 
real photonic decay width cannot be used to fix the coupling constants of 
these terms. The Dalitz decay width calculated from these terms is signifi- 
cantly smaller than the other contributions throughout the spectrum (with 
the exception of the 1/2— case at large dilepton masses), therefore they 
would need a much larger coupling constant to give a comparable contribu- 
tion. There seems to be no natural explanation for such a huge difference 
of the coupling constants. In the following we discard the terms discussed 
above and make the assumption that in the spin-1/2 case only the g\ term 
and in the spin>3/2 case only the gi and gi terms contribute to the Dalitz 
decay width. 
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Figure 1: Contributions of the various terms to the differential width of the 
Dalitz decay R — > Ne + e~ of an R resonance with mass = 1.5 GeV and 
with various spin-parities. The dimensionless coupling constants are set to 
unity. The explanation of the linestyles in the Figure for J p = 5/2— is valid 
for all plots for J >3/2. 
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Figure 2: Differential Dalitz decay width of a hypothetical baryon resonance 
with Breit-Wigner mass m^ w = 1.5 GeV, full width T to t = 0.12 GeV and 
a photonic branching ratio of 5%. The two linestyles correspond to the 
assumptions that only the g\ or g 2 terms contribute (g2 or gl ps 0). 
Predictions for resonances with masses above and below the Breit-Wigner 
mass (m* = mf w ± T tot ) are also shown. 
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Figure 3: Differential Dalitz decay width of hypothetical baryon resonances 
with mass m, = 1.5 GeV, photonic width T R ^ Nl = 0.006 GeV and with 
various spin-parities. 

In Fig. ^| the differential Dalitz decay width of hypothetical resonances 
with Breit-Wigner mass mf = 1.5 GeV, full width T to t = 0.12 GeV and a 
photonic branching ratio of 5% is shown for various spin-parities (spin>3/2). 
The two linestyles refer to the results obtained assuming that only the g\ or gi 
term contributes {gi ~ or g\ ^ 0). We also show the results for resonances 
with considerably larger and smaller mass, specifically, m* = mf w + T tot 
and m* = mf w — T tot . The real photonic width was used to fix the coupling 
constants in the m* = mf w case and the same coupling constants have been 
used in the m* = mf w ± T tot cases. 

For m* = mf w the gi and gi curves nearly coincide showing that these 
contributions are practically indistinguishable in the dilepton invariant mass 
spectrum if the mass of the decaying resonance is close to the Breit-Wigner 
mass. However, in the spin-parity 3/2+ and 5/2— cases the difference be- 
tween the gi and gi contributions is sizable if the resonance mass is well above 
the Breit-Wigner mass. Since in transport models masses of resonances are 
usually not far from the Breit-Wigner mass, transition from one of the g\ and 
gi terms to the other or to a linear combination will not have a substantial 
effect on the predictions for the dilepton invariant mass spectrum. In the 
following we make the assumption that only the g\ term contributes to the 
Dalitz decay width of baryon resonances. The explicite expressions for the 
polarization averaged squared transition matrix elements calculated from the 
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gi terms are given in Appendix B. 

To demonstrate the significance of spin and parity we show in Fig. ^ 
the differential Dalitz decay width of hypothetical resonances with mass m* 
= 1.5 GeV, photonic decay width = 0.006 GeV and with various 

spin-parities. The invariant mass spectra show large differences in the large 
dilepton mass range. The largest contribution comes from spin-parity 1/2+ 
and 3/2— resonances. 

The differential Dalitz decay width of real baryon resonances is shown 
in Fig. |j. We included all spin<5/2 nucleon and A resonances below 2 
GeV listed in the Review of Particle Physics [O] with at least *** status. 
We assumed that g 2 = g% = and fitted the g± coupling constants to the 
R + — > decay width, the numerical values of which we obtained using the 
Breit-Wigner full width and the pj branching ratio of the resonances given 
HI]. For the branching ratio, which is poorly known in many cases, we 
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used the average of the lower and upper bounds. In the case of A resonances, 
where the pj and n / y branching ratios are not listed separately, we used the 
A^7 partial width to fix the coupling constants. Note that the pj andn7 
branching ratios may differ significantly. 

In Table [I] we give the resulting coupling constants together with the 
values of the full width and the photonic branching ratio used in the fit. In 
the case of nucleon resonances we also give the couplings. 



4 Conclusion 

In this paper we studied the Dalitz decay of baryon resonances. We have 
given expressions for the differential Dalitz decay width of resonances with 
spin<5/2. The number of independent terms in the electromagnetic tran- 
sition matrix element is 2 in the spin-1/2 case and 3 in the spin>3/2 case. 
One of these terms does not contribute to the real photonic decay of the 
resonances, which is used to fix the coupling constants. We have shown that 
the remaining two matrix elements in the spin>3/2 case provide very similar 
contributions to the dilepton invariant mass spectrum. On the other hand, 
expressions for different spin-parities give very different dilepton spectra in 
the large dilepton mass region. Our results can be used in models of dielec- 
tron production in elementary hadronic reactions or heavy ion collisions. 
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Figure 4: Differential Dalitz decay width of unffavored baryon resonances. 
The curves correspond to resonances with m* = mf w . The Dalitz decay 
width of resonances with a mass different from the Breit-Wigner mass can 
be significantly different (c.f. Fig. 0). 
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J p 


resonance 


^tot 


Branching ratios (%) 


Coupling constants 

gV gV 


1/2+ 


N(1440) 
N(1710) 
A 1910) 


0.35 

0.1 

0.25 


0.042 0.021 
0.026 0.01 
0.1 


0.139 0.098 
0.030 0.019 
0.066 


1/2- 


N(1535) 
N(1650) 
A(1620) 


0.15 
0.15 
0.15 


0.2 0.15 
0.11 0.09 
0.024 


0.634 0.549 
0.315 0.285 
0.162 


3/2+ 


A(1232) 
N(1720) 
A(1600) 


0.12 
0.15 
0.35 


0.56 
0.05 0.2 
0.01 


1.98 

0.193 0.386 
0.162 


3/2- 


N(1520) 
N(1700) 
A(1700) 


0.12 

0.1 

0.3 


0.51 0.42 
0.03 0.07 
0.19 


0.793 0.719 
0.126 0.193 
0.549 


5/2+ 


N(1680) 
A(1905) 


0.13 
0.35 


0.27 0.034 
0.02 


2.74 0.971 
0.713 


5/2- 


N(1675) 
A(1930) 


0.15 
0.35 


0.014 0.07 
0.01 


0.678 1.52 
0.479 



Table 1: Coupling constants gi for the electromagnetic transitions of baryon 
resonances. In the case of nucleon resonances we give separately the coupling 
constants for the p7 and rry channels. We also give the full widths and 
photonic branching ratios used to obtain the values of the coupling constants. 
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Appendix A 

We give here the spin projectors needed to calculate the square of the tran- 



sition matrix elements. Throughout the paper we use the conventions of [12 



for relativistic quantities and Dirac matrices and spinors. Specifically, we use 
the normalization u(p, s)u(p, s) = 2m for Dirac spinors (p 2 = m 2 ), therefore 



y~]u(p, s)u(p, s) = f 1 + m. 

s 

The spin-3/2 projector: 

s)u u (p, s) = (p 1 + m) 

s 

The spin-5/2 projector: 

J2u^(p,s)u p °(p,s 

s 

3 



2 p^p u jAy" — p u Y^ 



3 m 2 



+ 



3m 



+ m) 
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{G w G va + G^G 



G^G pa 



with 



and 



2 17 7 7 7 J f 2m2 2m2 



Appendix B 
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(13) 



(14) 
(15) 

(16) 



In Eq. (^|) the squared transition matrix element averaged over initial and 
summed over final polarization states can be written as 



1 



n pol,R p i 



ql i q u 



9 + , 
Q 



M 



/IV 1 



(17) 
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where the notation 



M 



fJ,V 



Y: (N\J,\R)(R\J u \N) 



n pol,R pol NtR 

was used. Following || we introduce the notations 



(18) 



T, 



(19) 



M.^ v can be decomposed in terms of transverse and longitudinal parts, M.t 
and Ml, in the form 



M 



M T T. V - M L L 



(20) 



In terms of Air and M. l the polarization averaged squared transition matrix 
element is expressed as 



1 



n pol,R p i 



Y,\(N7\T\R)\ 2 = 2M T + M L . 



(21) 



Here we give the contributions of the g\ terms to M.t and M. l for the various 
spin-parity cases: 



M 1 ' 2 - 



Mll 2 - 



mV 2 ~ 



M 3 J 2+ 



Anag 2 



Arrag 1 
Anag 2 
Anag 2 



2m 4 
2 M_ 2 

2m 4 
1 



2 2 

— m 



(m, + mf - M 2 



ATragf-^-^im^+m) 2 (m* + m) 2 - M 2 



2m 4 
M 2 

2~^ 



(ml 



— m 



,m* - m 



(m, - m) 2 - M 2 
fm, - m) 2 - M 2 



(22) 
(23) 
(24) 
(25) 



12m 2 m 2 



(m* - m) 2 - M 2 



377^ + 6m„m + Amjn + 2m*m + 
-2m*mM 2 - 2m 2 M 2 + M 4 



(26) 
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M 3 / 2+ 



M% 2 - 



= 4nag 2 
= 4nag 2 



M b ' 2+ 



M] /2+ 



M 2 
3m 2 



(m* - mf - M 2 

m* + mf - M 2 



\2m 2 m 2 

(a 3 2 2 S A 

x \3m^ — 6m^m + Am^m — 2m*m + m 
+2m*mM 2 - 2m 2 M 2 + M 4 ) , 
M 2 r 

Airagf—r (m* + m) 2 — M 2 
3m A L 



(m* - m) 2 - M 5 



,772* + m) 
,4 



480m4m 4 

/a q *3 

x (2m* — Am^m + 3mjn - 2m»m + m 

+2m # mM 2 - 2m 2 M 2 + M 4 ) , 
/If 2 

m* — m) — M (to* + to) 



Anag 1 



47rag 1 



120to 2 to 4 
1 

'to* - m) 2 - M 2 



(27) 



(28) 
(29) 



M 



— M z 



(30) 
,(31) 



480mim 4 

x (2to 4 + Am\m + 3to 2 to 2 + 2m t m 3 + to 
-2m,mM 2 - 2m 2 M 2 + M 4 ) , 



+ m) 2 - ikP 

4 



= Ajtag 2 



M 2 



120m 2 m 4 



(m* — m) 2 — M 2 (m* + m) 



M 



(32) 
.(33) 
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